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ABSTRACT 


As part of the Seven Southeast Asian Studies (7SEAS) program, an Aerosol Robotic Network 
(AERONET) sun photometer and a Micro-Pulse Lidar Network (MPLNET) instrument have been 
deployed at Singapore to study the regional aerosol environment of the Maritime Continent (MC). 
In addition, the Navy Aerosol Analysis and Prediction System (NAAPS) is used to model aerosol 
transport over the region. From 24 September 2009 to 31 March 2011, the relationships between 
ground-, satellite- and model-based aerosol optical depth (AOD) and particulate matter with 
aerodynamic equivalent diameters less than 2.5 um (PM2>5) for air quality applications are 
investigated. When MPLNET-derived aerosol scale heights are applied to normalize AOD for 
comparison with surface PM2>s5 data, the empirical relationships are shown to improve with an 
increased 11 %, 10 % and 5 % in explained variances, for AERONET, MODIS and NAAPS 
respectively. The ratios of root mean square errors to standard deviations for the relationships also 
show corresponding improvements of 8 %, 6 % and 2 %. Aerosol scale heights are observed to be 
bimodal with a mode below and another above the strongly-capped/deep near-surface layer (SCD; 
QO — 1.35 km). Aerosol extinctions within SCD are well-correlated with surface PM25 
concentrations, possibly due to strong vertical mixing in the region. 


Keywords: Air Pollution, Air Quality, Aerosol Optical Depth 
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1 INTRODUCTION 


In the form of respirable particulate matter with aerodynamic equivalent diameter of less 
than 10 um (PMio) and 2.5 um (PM2:5), atmospheric aerosol particles pose a significant public 
health risk. Common health implications include cardiopulmonary conditions, acute and chronic 
respiratory infections, and lung cancer (Cohen et al., 2005). In response to worsening regional air 
pollution in Southeast Asia (SEA) and the Maritime Continent (MC), due to increased industrial 
and urban emissions (e.g., Salinas et al., 2009; Reid et al., 2013), as well as regional biomass 
burning smoke (e.g., See et al., 2006, 2007; Salinas et al., 2013a, 2013b) from land preparation and 
forest clearing activities (e.g., Field et al., 2009), several ground-based air quality monitoring 
networks, such as those operated by Malaysia’s Department of Environment (DOE; 
http://apims.doe.gov.my), Singapore’s National Environment Agency (NEA; 
http://www.nea.gov.sg) and  Thailand’s Pollution Control Department (PCD; 
http://www.pcd.go.th), have been established. These networks provide information on air quality, 
including mass concentrations for PM, carbon monoxide, nitrogen oxide, sulphur dioxide and 
ozone, which serves as the basis for issuing public health advisories in cases of deteriorating air 
quality. While such networks are suitable for assessing air pollution from local sources, they are 
relatively limited spatially, and have poor representativeness of the vast rural areas in the region. 
Therefore, satellite-based remote sensing observations (e.g. Hoff and Christopher, 2009) and 
aerosol transport models (e.g., Hyer and Chew, 2010; Wang et al., 2013; Kim et al., 2014) with 
wide spatial coverage can be used to complement the operations of air quality monitoring networks, 
and provide needed information on the spatial distribution of regional air pollution, especially in 
cases of significant biomass burning outbreak (e.g., Reid et al., 2013). 

Chu et al. (2003), as well as Wang and Christopher (2003), first demonstrated the 
capabilities of the Moderate Resolution Imaging Spectroradiometer (MODIS; Remer et al., 2005) 
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on board the National Aeronautics and Space Administration’s (NASA) Terra and Aqua satellites 
to correlate satellite-based aerosol optical depth (AOD) with surface PM19 and PM2.5 measurements 
respectively. Liu et al. (2004; 2005) also showed the Multi-Angle Imaging Spectroradiometer 
(MISR; Diner et al., 2002) on board Terra to have similar capabilities for air quality applications. 
Several research efforts have worked towards establishing empirical satellite-based AOD-PM 
relationships in different regions of the world (e.g., Engel-Cox et al., 2004b; Gupta et al., 2006; 
Schaap et al., 2008), as well as deriving satellite-based AOD products with higher resolution for 
comparison to urban domains (e.g., Li et al., 2005; Kumar et al., 2007; Wong et al., 2011; Bilal et 
al., 2013). Meteorological parameters, such as relative humidity (RH), wind direction and speed 
(e.g., Gupta and Christopher, 2009a, 2009b), are also shown to be important in establishing AOD- 
PM correlations. 

There are certainly large-scale relationships between ambient columnar integrated AOD 
and surface PM concentration, but beyond this, the variability in aerosol vertical distribution and 
particle hygroscopicity lead to significant regional variability between the two parameters (Toth et 
al., 2014). Although Engel-Cox et al. (2006) has qualitatively demonstrated the potential of lidar- 
derived aerosol vertical distribution to constrain the AOD-PM relationship, the increasing 
availability of continuous long-term global ground-based lidar measurements provides the 
necessary data to assess the quantitative impact (Tsai et al., 2011; Chu et al., 2013). Given the 
complex meteorological environment of SEA and the MC, such an analysis requires the 
incorporation of continuous long-term lidar measurements over many seasons to limit inter-annual 
variance (Campbell et al., 2013; Reid et al., 2012, 2013). 

Singapore Aerosol Robotic Network (AERONET; Holben et al., 1998) and Micro-Pulse 
Lidar Network (MPLNET; Welton et al., 2001) datasets have previously been considered in tandem 


for quantifying and characterizing the seasonal variability of aerosol particle columnar optical 
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properties and vertical distributions at the National University of Singapore (NUS) atmospheric 
measurement supersite from 24 September 2009 to 31 March 2011 (Chew et al., 2013). In this 
study, modified aerosol scale heights (Hayasaka et al., 2007) are derived from vertically-resolved 
MPLNET AOD profiles, and used to evaluate AERONET, MODIS and Navy Aerosol Analysis 
and Prediction System (NAAPS; Witek et al., 2007; Reid et al., 2009; Westphal et al., 2009) 
modeled AOD-PM2.s5 relationships over Singapore for the same study period. Over Singapore, 
transported aerosol layers within the residual boundary layer may be mixed into the surface layer 
by way of the entrainment zone (Atwood et al., 2013). This mechanism might complicate the 
development of any AOD-PM 25s (empirical or otherwise) correlations. The correlations will also 
likely be affected by varying aerosol hygroscopic growth coupled with high RH, resulting in yet 
more variable light scattering efficiencies in ambient conditions (Salinas et al., 2013a; 2013b). 

In this paper, we describe a study to investigate the extent to which satellite, ground based 
and model simulated AOD relates surface PM in Singapore from 24 September 2009 to 31 March 
2011. This is conducted with the aid of the aforementioned MPLNET instrument to help interpret 
aerosol vertical distributions. The modified aerosol scale heights (e.g., Hayasaka et al., 2007) is 
derived from MPLNET data, and used to normalize AOD in order to estimate the hourly surface 
PM2s5 concentrations. A simple hygroscopic growth correction based on the assumption of a sulfate 
environment in Singapore (See et al., 2006) is also applied to minimize the differences between 
ambient PM2.s5 concentrations and the measured PM2.5 dry mass. Finally, in order to evaluate the 
models of PM2.5 estimation, the improvements in explained variances, the root-mean-square error 


(RMSE), normalized RMSE (NRMSE) and ratio of RMSE to standard deviation (RSR) are used. 


2 INSTRUMENTS, DATA AND METHODS 


2.1 PM)2.s Measurements 
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NEA operates thirteen air quality monitoring stations in Singapore (NEA, 201 1a), of which 
Pandan Reservoir station is nearest to the NUS atmospheric measurement supersite (~ 4 km; For 
supersite details, see Atwood et al., 2013; Chew et al., 2013; Salinas et al., 2013a). These stations 
are positioned strategically in industrial, urban and suburban areas, as well as along roadsides to 
monitor the ambient air quality. They are linked up via a telemetric monitoring and management 
system. The concentrations of major criteria pollutants (i.e. PM25, PMio, carbon monoxide, 
nitrogen oxide, carbon monoxide, and sulphur dioxide) are measured, and reported to the public 
daily. For this study, PM2.5 measurements are applied as hourly averages. 

At Pandan Reservoir, PM2.5 concentrations are measured with the Thermo Scientific FH 62 
C14 series continuous ambient particulate monitor. Air is first sampled via an aerosol inlet with 
cutoff diameter of 2.5 um, which is heated above ambient temperature to remove condensation. 
Sampled particles are subsequently collected onto a filter tape between a Carbon-14 source and 
detector within the instrument. The instrument measures the concentration of sampled particles via 
the beta ray attenuation method (e.g., Hinds, 1999). The filter tape sample spot is advanced every 
hour, thus providing hourly PM2.5 mass concentrations. 

Due to the high RH in the region, aerosol hygroscopic growth influences the AOD-PM25 
relationship (Tsai et al., 2011; Chu et al., 2013). When particles containing hygroscopic 
components (e.g., sulfates and some organics) are exposed to high RH (> 70%), the scattering 
cross-section of particles undergoing substantial water absorption can increase by more than a 
factor of two (e.g., Hanel 1976). Even biomass burning particles, which are globally much less 
hygroscopic than anthropogenic emission counterparts, have a high hygroscopicity due to their 
high primary and/or secondary sulphur load — a result perhaps of the regions rich in volcanic soils 
(e.g., Reid et al., 2005; 2013). Based on the physics of highly soluble salts, the hygroscopic growth 


function is frequently written as 
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ARH) = (RA oo) ° 


(1) 


where RH is in percentage, RHo is referenced at 30 % and g is an empirical growth fitting parameter 
(Kasten, 1969). Values of g range from zero for insoluble particles to almost unity for highly 
soluble particles. 

Due to the lack of aerosol hygroscopicity measurements from a dual nephelometer set-up 
(e.g., Hegg et al., 1996; Kotchenruther and Hobbs, 1998), which renders more constrained 
estimates, the hygroscopic growth function of a sulfate environment is assumed (i.e., g = 0.63; 
Hanel, 1976) and applied to the PM2.5 measurements here. This assumption is supported by the 
observations of See et al. (2006) in reporting in situ measurements of aerosol microphysical 
properties at Singapore as consistent with that for a sulfate pollution environment (e.g., Waggoner 
and Weiss, 1980). 

Hourly-averaged RH data are acquired from the National Climate Data Center (NCDC; 
http://www.ncdc.noa.gov), and used to compute f(RH). Although the annual average RH is 84.2%, 
Singapore experiences a diurnal variation from ~ 90% in the early morning to 60% in the mid- 
afternoon, even going below 50% at times (NEA, 2011b). During prolonged heavy rain, RH often 
approaches 100%. As substantial errors in the calculation of aerosol hygroscopic growth factor can 
be introduced when aerosol types are not properly identified, only data with RH < 75% are used 


(e.g. Tsai et al., 2011; Chu et al., 2013) in this study. 


2.2. AERONET and MPLNET 
The Singapore AERONET site deploys a Cimel Electronique CE-318N sun photometer. 


The instrument collects direct solar measurements in eight spectral bands (with center wavelengths 
6 


160 


161 


162 


163 


164 


165 


166 


167 


168 


169 


170 


171 


172 


173 


174 


175 


176 


1a 


178 


179 


180 


181 


182 


from 0.340 um to 1.640 um) every 30 s within one-minute periods, which are then averaged as 
triplet measurements in order to derive AOD (Holben et al., 1998). Only cloud-screened (Smirnov 
et al., 2000) and quality-assured Level 2.0 AOD data are used in this study. The estimated accuracy 
of AERONET Cimels varies spectrally from +0.01 to +0.02 for AOD measurements with higher 
errors in the ultraviolet channels (Holben et al., 1998; Eck et al., 1999). Chew et al. (2011) found, 
however, that unscreened cirrus clouds significantly increase this relative uncertainty at Singapore, 
with potential cirrus artifact of 0.02 to 0.08 being common. 

The collated Singapore MPLNET instrument (0.527 um; Spinhirne, 1993; Spinhirne et al., 
1995) is acompact and eye-safe mie lidar. The instrument is able to profile aerosol and cloud layers 
vertically by transmitting short laser pulses into the atmosphere and then measuring the 
backscattered signals. Instrument calibration includes corrections for detector dead-time, dark 
counts, afterpulsing and lidar overlap (Campbell et al., 2002). After subtracting the solar 
background signal, the raw profiles are then converted to MPLNET Level 1.0 Normalized Relative 
Backscatter (NRB) product, which is range- and energy-normalized. The uncertainties for the NRB 
product have been detailed in Welton and Campbell (2002). 

MPLNET Level 1.0 NRB product are first screened for visible clouds (e.g., Clothiaux et 
al., 1998) and optically-thin cirrus clouds (e.g., Campbell et al., 2008; Chew et al., 2011; Campbell 
et al., 2015b), and subsequently averaged as 20-minute segments about AERONET Level 2.0 AOD 
observations. The former step helps filter potential error in collocated AERONET observations 
from unscreened cloudiness. Vertical profiles of aerosol particle backscatter and extinction 
coefficients are derived from these 20-minute averaged NRB data with the backward Fernald two- 
component solution to the lidar equation (Fernald, 1984; Welton et al., 2000). These profiles are 


considered quality-assured as MPLNET Level 2.0a aerosol product when they are acquired within 
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an optimal temperature range (23 + 5 °C for the Singapore MPLNET instrument), have extinction- 
to-backscatter ratio error < 30%, and include > 80% of 20-minute averaged NRB data. 

In this study, modified aerosol scale heights (/7,,) are derived from the MPLNET Level 2.0a 
vertical profiles of extinction coefficient (f(z) as a function of height z) such that approximately 


63% of the total AOD exist below the point (e.g., Hayasaka et al., 2007; Turner et al., 2001): 


Hm 
B(z) dz =(1—e7-1)AOD = 0.632 AOD. (2) 


0 


H,, is subsequently used to normalize AOD in order to estimate surface PM25 concentrations. 
Although mean planetary boundary layer (PBL) extinctions can also be used to estimate surface 
PM2.s concentrations (e.g., Engel-Cox et al., 2006; Liu et al., 2011), the derivation of PBL heights 
over the region can be erroneous due to the frequent presence of low-level clouds causing deeper 
PBL height retrievals with the typical wavelet techniques (e.g., Brooks, 2003). At this point of 
time, MPLNET is processing all data, including that of Singapore, which will be collectively 
known as “Version 3” and includes an improved PBL height product (Welton, 2015; personal 
communication). The new PBL height retrieval algorithm uses image processing and fuzzy logic 
techniques to overcome the shortcomings of using wavelet techniques alone (Lewis et al., 2013). 
Fig. 1 depicts seasonal H,,, distributions over Singapore from October to December 2009 
(OND 2009; Fig. la), January to March 2010 (JFM 2010; Fig. 1b), April to June (AMJ 2010; Fig. 
1c), July to September (JAS 2010; Fig. 1d), October to December 2010 (OND 2010; Fig. le) and 
January to March 2011 (JFM; Fig. 1f), corresponding with the seasonal aerosol extinction 
coefficient profiles, lidar ratios and Angstrom exponents reported in Chew et al. (2013). A bi-modal 


distribution with a mode below and above the top of strongly-capped/deep near-surface layer (SCD; 
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0 — 1.35 km; Chew et al., 2013) is found, consistent with the presence of vertical wind shear over 
Singapore, resulting in potential aerosol reservoirs aloft above the boundary layer (Atwood et al., 
2013). The overall H,,, distribution exhibits a mean of ~ 1.4 km, as summarized in Fig. 1g, with 
42% of occurrence below and 58% of occurrence above the top of the SCD. 

Yuet al. (2010) found CALIOP-derived seasonal H,,, of ~ 1.4 — 1.9 km for the SEA region, 
generating seasonal H,,, of ~ 1.7 —2.1 km with the Goddard Chemistry Aerosol Radiation Transport 
(GOCART) model. While CALIPO-derived H,, are more consistent with the MPLNET-derived 
quantities, GOCART-derived quantities are consistently higher than ground- and satellite-based 


measurements, indicting a likely overestimation of aerosol transport within the free atmosphere. 


2.3. MODIS Level 2 Aerosol Products 

MODIS Collection 5.2 aerosol products (Level 2 MOD04 and MYDO4 for data collected 
from NASA Terra and Aqua satellite respectively) are available from the NASA Level | and 
Atmosphere Archive and Distribution System (LAADS). Alternatively, the data can be generated 
with the International MODIS/AIRS Processing Package (IMAPP; 
http://cimss.ssec.wisc.edu/imapp/) for ground stations capable of receiving direct broadcast from 
the Terra and Aqua satellites, which can be advantageous for monitoring air pollution in near real 
time. For this study, the AOD measurements at 0.55 um (Optical_Depth_Land_And_Ocean 
variable of MOD04 and MYDO4 aerosol product) with a nominal spatial resolution of 10 by 10 
km (Remer et al., 2005) are used. The actual ground footprint for these satellite AOD data ranges 


from 10 by 10 km at nadir to roughly 40 by 20 km at the edge of the satellite swath. 


2.4 NAAPS aerosol transport model and FLAMBE smoke flux model 
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The NAAPS model is a modified form of a hemispheric model of sulfate aerosols developed 
by Christensen (1997), with dust, biomass burning smoke and sea salt as described by Westphal et 
al., 2009; Reid et al., 2009; and Witek et al., 2007, respectively. For this study period, NAAPS is 
driven by meteorological fields from the Navy Operational Global Atmospheric Prediction System 
(NOGAPS; e.g., Hogan and Rosmond, 1991; Hogan and Brody, 1993). Smoke, dust and sulfate 
AOD, as well as aerosol mass concentrations, are determined at a 1° x 1° grid for this study in 6-h 
intervals and twenty-four vertical levels up to 100 mb. The Naval Research Laboratory (NRL) 
Atmospheric Variational Data Assimilation System — Aerosol Optical Depth (NAVDAS-AOD; 
Zhang et al, 2008) is used to improve NAAPS model simulations by assimilating MODIS AOD 
data over sea and land within the modeled fields (Hyer et al., 2011; Shi et al., 2011). 

The Fire Locating and Monitoring of Burning Emissions (FLAMBE) smoke flux model is 
applied to estimate biomass burning emissions within NAAPS (Reid et al., 2004; 2009). MODIS 
active fire detections are used to determine the locations and timings of fires (Giglio et al., 2003), 
which are subsequently mapped onto the Global Land Cover Classification (GLCC Version 2.0; 
Loveland and Belward, 1997) 1-km land cover map to estimate burning fuels. The resultant fire 
intensive properties (i.e., fuel loading, fuel consumption, emission partitioning according to various 
aerosol species) are then calculated based on the fuel models of Reid et al. (2005). 

Hyer and Chew (2010) have previously used NAAPS analyses to examine an extreme MC 
biomass burning smoke episode from September to October 2006. In that study, NAAPS was able 
to capture the timing of smoke events, but severely underestimated their magnitude. Campbell et 
al. (2015) recently conclude a similar result, having paired FLAMBE with WRF-Chem and 
examining 3-5 day smoke transport across Borneo Island. These studies indicate that FLAMBE is 
not fully resolving the variability in smoke emissions. Though, based on MODIS overpasses alone, 


this is to be expected given only four possible views of a given scene per day and only a coarse 
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mapping of available fuel sources. Based on a regression analysis between the measured daily PM10 
concentrations and NAAPS surface concentrations at 54 air quality monitoring stations in 
Singapore and Malaysia, we estimate that FLAMBE is underestimating regional smoke emissions 
by a factor of between 2.5 to 10, and a overall correction factor of 3.5 is approximated. This 
correction factor is similar to other regional studies conducted at the mesoscale (e.g., Wang et al., 


2013). 


3 RESULTS AND DISCUSSION 


3.1 Correlation between MPLNET-derived extinction coefficients and hygroscopic- 

corrected PM2.5 measurements 

After constraining our data sample for cases of RH < 75%, there are a total of 855 
AERONET Level 2.0 data coincident with cloud-free MPLNET-derived vertical profiles of aerosol 
particle extinction available for analysis. In order to match the temporal resolution of the PM2.5 
measurements, available AERONET and MPLNET data are also hourly averaged (N = 205). 

MPLNET-derived extinction coefficients at 100-m intervals from 0 to 2.4 km are compared 
with hygroscopic-corrected surface PM2.5 concentrations, with the resultant coefficients of 
determination (R’) shown in Table 1.It should be noted that overlap corrections can be significant 
for heights lower than 200 m (Campbell et al., 2002; Toth et al., 2014). Correlations between the 
two are moderate (R” ~ 0.35 — 0.54) from 0 — 1.2 km, which includes the SCD. Above 1.3 km, 
correlation between the two is substantially reduced (R* ~ 0.03 — 0.29; not shown). Higher 
correlation within the SCD may be the result of strong vertical mixing produced by updrafts due to 
high surface temperatures over Singapore (e.g., Li et al., 2013). As such, a mean extinction 


coefficient within the mixed layer is likely more suitable than a series of contiguous near-surface 
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extinction coefficients for comparison with hourly-average PM2s5 measurements, taking into 


account possible vertical mixing within an hourly time frame. 


32 Correlation between AERONET AOD and _=shygroscopic-corrected PM2.s5 
measurements 

H,,, derived from MPLNET Level 2.0a data are used to scale AERONET, MODIS and 
NAAPS AOD to approximate a near surface extinction coefficient for comparison with ground- 
level PM2.s measurements, similar to the methods described in Tsai et al. (2011) and Chu et al. 
(2013). Substantial variability in aerosol vertical distributions exists over Singapore (Chew et al., 
2013), due primarily to vertical wind shear and the potential for aerosol reservoirs aloft above the 
boundary layer (e.g., Atwood et al., 2013). Therefore, scaling the columnar AOD by the modified 
aerosol scale height is a simple method that takes into account the possible air mass exchange 
between the boundary layer and free troposphere. 

Fig. 2 features correlation between AERONET AOD and hygroscopic-corrected PM2.5 
concentration as aerosol scale height corrections are applied. Correlation is high (R? ~ 0.47 — 0.58), 
with an increase of 11% in explained variances (from the difference in R* values) when aerosol 
scale height corrections are applied. Scatter about the regression line is also noticeably reduced. 

In order to further evaluate the models of PM25 estimation by AOD and AOD/H,, , we 
calculate RMSE, NRMSE and RSR for the AERONET data in Table 2, together with MODIS and 
NAAPS data for comparison (e.g., Chu et al., 2013). RMSE indicates the difference between the 
fitted and observed data, whereas NRMSE also considers the range, i.e. the maximum and 
minimum, of the observed data. RSR further compares the uncertainty of PM2.s5 estimated from 
AOD and AOD/H,,, to the observations. RMSE (10.32 — 11.61 wg m™) is comparable with NAAPS 
data, but NRMSE (0.10 — 0.11) is the lowest among all the data. The hygroscopic-corrected PM2.5 
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as estimated by AERONET data has the lowest RSR (0.65 — 0.73) among all the data, as well as 


an improvement of 8 % in RSR when using scaled AERONET AOD. 


3.3 Correlation between MODIS AOD and hygroscopic-corrected PM2.5 measurements 

In order to spatially correlate ground measurements (i.e. MPLNET-derived H,,, and PM2:5 
concentrations) with satellite-based MODIS AOD, they are paired up with a satellite observation 
by choosing image pixels that are spatially and temporally close. First, only ground measurements 
within an hour of the satellite observations are considered for correlation. Second, all satellite 
observations within 0.50° (~ 40 — 50 km) of the measurement supersite are spatially averaged, due 
to the constant cloud coverage over Singapore preventing otherwise routine data acquisition. Such 
separation requirements are consistent with the lagged correlation analysis conducted for surface 
measurements against sky-based and space-based measurements by Anderson et al. (2003). 

Although 220 MODIS AOD observations are available from September 2009 to March 
2011, only 38 observations are available for comparison with ground measurements due to 
separation requirements. Fig. 3a illustrates the relatively tenuous relationship between AERONET 
and MODIS AOD (R? = 0.54). Despite small differences in sampling methodology, Hyer et al. 
(2011) found varying slopes from 0.76 (Mukdahan; 17° N, 105° E) to 1.52 (Bac Lieu; 9° N and 
106° E) when comparing between AERONET and MODIS AOD in SEA due to the complex 
aerosol and land surface environment. 

Fig. 3b and c depict correlation between MODIS AOD and hygroscopic-corrected PM2:s 
concentration as aerosol scale height corrections are applied. The correlation for MODIS is weaker 
when compared with AERONET, although improvements (from R? of 0.35 to 0.45) are noted with 
an increase of 10% in explained variances (from the difference in R* values) when aerosol scale 


height corrections are applied. The primary reason for the lesser correlation is attributed to constant 
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cloud cover and contamination over the region, resulting in less quality-assured MODIS AOD 
retrievals available near Singapore. Spatial averaging of MODIS AOD to compensate for cloud 
cover is also likely to reduce representativeness. 

RMSE, NRMSE and RSR are also calculated for MODIS data in Table 2. Although RMSE 
(7.90 — 8.55 ug m*®) is the lowest among all the data, this result is likely an artifact due to the small 
sample size. NRMSE (0.22 — 0.23) for MODIS data is the largest among all the data, and is 
indicative that RMSE is large compared to the observed data range. When using MODIS AOD/H,, 
instead of AOD to estimate hygroscopic-corrected PM2.s, there is an improvement of 6 % in RSR 


(from 0.82 to 0.76). 


3.4 Correlation between NAAPS AOD, surface concentrations and hygroscopic-corrected 

PM2.5 measurements 

As satellite observations over Singapore are limited by relatively frequent cloud cover, 
NAAPS AOD can be considered a first-order proxy for MODIS AOD, and has in fact been shown 
to correlate well with MODIS in Southeast Asia in the 00 hr model analysis state through 
verification with AERONET (Campbell et al., 2013). In order to spatially correlate NAAPS AOD 
and surface concentration with ground measurements (i.e. MPLNET-derived H,, and PM25 
concentrations), only output from the nearest NAAPS model grid are used and interpolated with 
regards to time. There exists an extreme outlier when comparing the AERONET and NAAPS AOD 
(not shown), which is subsequently removed. 

Fig. 4a features correlation between NAAPS surface concentration and PM25 
concentration. Despite AOD data assimilation, NAAPS underestimates surface level 
concentrations over Singapore, as previously observed by Hyer and Chew (2010) for PMio 


concentrations. NAAPS representation of near surface processes is limited by the coarse model 
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spatial resolution and rapid vertical mixing scheme among the lowest model layers, thus aerosol 
vertical distributions and surface concentrations are not realistically simulated. 

In contrast, AOD is a columnar optical quantity, which is not affected by aforementioned 
model limitations. Fig. 4b and c reflect reasonable correlation between NAAPS AOD and 
hygroscopic-corrected PM2s concentration (R* ~ 0.22 — 0.27) with and without aerosol scale height 
corrections, respectively. However, any improvement in correlation from applying the aerosol scale 
height corrections is considered marginal, as there is only 5% increase in explained variance (see 
difference in R? values). 

RMSE, NRMSE and RSR are calculated for NAAPS data in Table 2. RMSE (12.82 — 13.24 
ug m3) is comparable to AERONET data. NRMSE (0.16) is similar to AERONET data, and lower 
than MODIS data due to the large sample size. There is also slight improvement of 2 % in RSR 
(from 0.88 to 0.86) when using scaled NAAPS AOD. These results are possibly due to model- 
related issues as discussed in Hyer and Chew (2010), especially the deficiencies in the emission 
database. A recent field experiment in Southeast Asia investigated MPLNET aerosol extinction 
profiles versus those derived from the Weather Research and Forecasting model coupled with 
Chemistry (WRF-Chem) using FLAMBE as the emissions (Wang et al., 2013; Campbell et al., 
2015a). Despite the higher model resolution allowed in WRF-Chem compared with NAAPS, 
model-derived aerosol extinctions are still mostly underestimated, as described above. Similarly, 
regional cloud cover prevents effective MODIS data acquisition and assimilation for NAAPS. 
Cloud bias in MODIS data can also be so severe in the region that it decreases model efficacy when 


unscreened observations are assimilated (Reid et al., 2009; Hyer et al., 2011). 


4 CONCLUSIONS 
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In the study, the relationships between Aerosol Robotic Network (AERONET), Moderate 
Resolution Imaging Spectroradiometer (MODIS), Navy Aerosol Analysis and Prediction System 
(NAAPS) aerosol optical depth (AOD) and particulate matter with aerodynamic equivalent 
diameter of less than 2.5 um (PM2>5) concentrations over Singapore for potential air quality 
applications are examined using coincident AERONET Level 2.0 and MPLNET Level 2.0a 
datasets from 24 September 2009 to 31 March 2011. The modified aerosol scale height (4,,,), which 
represents the vertical extent of aerosols (Hayasaka et al., 2007), are derived to scale columnar 
AOD to approximate near surface extinction for comparison with surface PM2s5 measurements. 
This approach is similar to the methods described in Chu et al. (2013) and Tsai et al. (2011), which 
normalized columnar MODIS AOD with a combined planetary boundary layer height (PBL) and 
a scale height of exponential decay of aerosol extinction above the PBL. The conclusions are as 


follow: 


1. The H,, distribution over Singapore is observed to be bimodal with a mode below and 
another mode above the top of strongly-capped/deep near-surface layer (SCD; 0 — 1.35 
km; Chew et al., 2013). The mean H,, of 1.4 km is derived, which is comparable to 
CALIOP-derived seasonal H,, of ~ 1.4 km — 1.9 km over the SEA region (Yu et al., 
2010). 

2. MPLNET-derived near surface extinctions (0 — 1.3 km) are well correlated with surface 
PMs concentrations (R? ~ 0.29 — 0.54) within the SCD layer, which may be due to 
strong vertical mixing caused by high surface temperature (Li et al., 2013). Above 1.3 
km, the correlations are substantially reduced (R ~ 0.03 — 0.21). 

3. The AERONET and MODIS AOD-PM_:; correlations are improved (i.e. from R* of 
0.47 to 0.58, and from R? of 0.35 to 0.45, respectively) when applying aerosol scale 
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height corrections, with an increase of 11% and 10% of explained variance respectively. 
There are also improvements of 8 % and 6 % when considering the ratio of root mean 
square error to standard deviation (RSR) for AERONET and MODIS AOD 
respectively. 

4. The NAAPS AOD-PM2 correlation changes marginally (i.e. from R? of 0.22 to 0.27) 
when applying aerosol scale height corrections, with a slight increase of 5% in 
explained variance. There is also a slight improvement of 2 % when considering RSR 
for NAAPS AOD. NAAPS surface concentration is also not well correlated with PM2.5 
concentrations. The weaker correlations are possibly due to the deficiencies in the 


FLAMBE database for the region. 


With the imminent launch of earth observation satellites at an near-equatorial orbit, such as 
Singapore’s Temasek Low Earth Orbit Satellite (TeLEOS-1) and Kent Ridge satellite (KR-1), as 
well as geostationary satellites, such as Himawari-8/9 
(http://www.data.jma.go.jp/mscweb/en/himawari89/index.html), more frequent passes over the 
region will provide timely monitoring of air pollution and observations of extreme biomass burning 
smoke. This will increase regional coverage, compared with sun-synchronous instruments like 
MODIS, and thus provide additional opportunity to alleviate data-denial conditions induced by 
endemic cloud cover. Further, greater satellite coverage, and corresponding aerosol observations, 
will allow better estimations of surface aerosol concentration, providing more data for civil 


authorities and more assimilable data for model systems. 
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FIGURE CAPTIONS: 

Figure 1. Seasonal modified aerosol scale height distributions detected by MPLNET instrument 
at 0.527 um over Singapore for (a) OND 2009, (b) JFM, (c) AMJ, (d) JAS, (e) OND 2010 and (f) 
JFM 2011. (g) Overall modified aerosol scale height distribution. Bin sizes are 0.2 km. Blue 


dashed line marks top of the SCD layer (0 — 1.35 km). 


Figure 2. Correlation between (a) AERONET AOD and _ hygroscopic-corrected PM2.5 
concentration, and (b) AOD/H,,, and hygroscopic-corrected PM2.s5 concentration. The text within 
each panel describes the coefficient of determination (R7), significance level of correlation (P, not 
significant when P greater than 0.05), number of comparison pairs (N), root-mean-square error 


(RMSE), and mean + standard deviation of measured quantities. 


Figure 3. (a) Validation between AERONET and MODIS AOD over Singapore. Correlation 
between (b) MODIS AOD and hygroscopic-corrected PM25 concentration, and (c) AOD/H,, and 
hygroscopic-corrected PM2.5 concentration. Text description within (b) and (c) is the same as for 


Figure 2. 


Figure 4. Correlation between (a) NAAPS surface concentration and PM2.5 concentration, (b) 


NAAPS AOD and hygroscopic-corrected PM2.5 concentration, and (c) AOD/H,, and hygroscopic- 


corrected PM2.s5 concentration. Text description within panels is the same as for Figure 2. 
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660 Table 1: Coefficients of determination (R?) between hygroscopic-corrected PM2s5 concentrations and MPLNET derived extinction coefficients at 
661  100-m intervals. It should be noted that overlap corrections can be significant at heights lower than 200 m. 
662 


100m 200 m 300 m 400 m 500 m 600 m 700 m 800 m 900m 1000m 1100m 1200m 


0.54 0.54 0.54 0.54 O55 0.57 0.53 0.46 0.38 0.40 0.46 0.35 


663 


664 Table 2: Number (N), mean, median, maximum (max) / minimum (min), and standard deviation (S.D.) for hygroscopic-corrected PM2.5 
665 measurements. Root-mean-square error (RMSE), normalized RMSE (NRMSB), and ratio of RMSE to standard deviation of observed quantities 


666 (RSR) for hygroscopic-corrected PM25 estimated by AOD and AOD/H,,. 


667 
nN 205 38 204 
on 24.33 2271 23.94 
HRESENE een 19.86 20.92 19.83 
Mini Max 1.53 / 103.00 6.97 / 43.68 1.53 / 82.80 
ae 15.92 10.46 14.97 
AERONET MODIS NAAPS 
er AOD AOD/H,, AOD AOD/H, AOD AOD/H,, 
einen RMSE 11.61 10.32 8.55 7.90 13.24 12.82 
NEMeE 0.11 0.10 0.23 0.22 0.16 0.16 
ROR 0.73 0.65 0.82 0.76 0.88 0.86 
668 
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